Dynamics of radiation due to vacuum nonlinearities 
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In quantum electrodynamics, photon-photon scattering can be the result of the exchange of vir- 
tual electron-positron pairs. This gives rise to a non-trivial dispersion relation for a single photon 
moving on a background of electromagnetic fields. Knowledge of the dispersion relation can be 
transferred, using standard methods, into new insights in the dynamical equations for the pho- 
tons. Effectively, those equations will contain different types of self-interaction terms, depending on 
~ ' whether the photons are coherent or not. It is shown that coherent photons are governed by a non- 

| linear Schrodinger type equation, such that the self-interaction terms vanish in the limit of parallel 

f») . propagating waves. For incoherent photons, a set of fluid equations can determine the evolution of 

■ the corresponding radiation gas. In the case of a self-interacting radiation fluid, it is shown that 
| t ' Landau damping can occur. 

'■ 

Q . PACS numbers: 12.20.Ds, 95.30.Cq 

t> : 

t-H ' I. INTRODUCTION 

> : 

According to QED, the non-classical phenomenon of photon-photon scattering can take place due to the exchange of 

■ virtual electron-positron pairs. This is a second order effect (in terms of the fine structure constant a = e 2 / Aireohc « 
1/137), which in standard notations can be deduced from the Euler-Heisenberg Lagrangian density 0,0 

m ■ JSf - £o^ + £o«( 4 ^ 2 + 7^ 2 ) , (1) 

O ; 

where n = 2a 2 h 3 / 45m^c 5 , & = \{E 2 — c 2 B 2 ), = cE • B, e is the electron charge, c the velocity of light, 2-k% the 
Planck constant and m e the electron mass. The Lagrangian Q is valid as long as there is no pair creation and the 

1 | field strength is smaller than the critical field, i.e. 



Oh 



X 



ui <C m e c 2 /h and \E\ <C -E C rit = m e c 2 /e\ c (2) 



respectively. Here A c is the Compton wave length, and E CT it ~ 10 18 V/m. The latter terms in JQl represent the effects 
of vacuum polarisation and magnetisation. We note that & = & = in the limit of parallel propagating waves. It is 
therefore necessary to use other waves in order to obtain an effect from the QED corrections. Several attempts have 
been presented in the literature H SHIES HQ, where Refs ' 0>SS 

6] mainly focused on principal issues, whereas 
the experimental possibilities for detection have been discussed in Refs. 0,SI3- 

The non-trivial propagation of photons in strong background electro magn etic fields, due to effects of nonlinear 
electrodynamics, has been considered in a number of papers (see, e.g., Refs. [lllll2j and references therein). Their main 
focus was on the interesting effects of photon splitting and birefringence in vacuum. However, Thoma [Lj investigated 
the interaction of photons with a photon gas, using the real-time formalism, and calculated the corresponding change 
in the speed of light due to the Cosmic Microwave Background (CMB). 

In the present paper, we derive an evolution equation for an ensemble of electromagnetic pulses. Moreover, the 
methods of radiation hydrodynamics 0] are combined with the QED theory for photon-photon scattering, and a 
system of equations (c.f. [l5l|) is obtained, where the radiation pressure of the pulse will act as a driver of acoustic 
waves in the photon gas. 
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II. EQUATIONS AND EXAMPLES 



In a medium with polarisation P and magnetisation M the general wave equations for E and B arc 



1 d 2 E 



V 2 E 



-Mo 



d 2 P 



■ c 2 V(V • P) 



d_ 

dt 



(V x M) 



(3a) 



and 



1 d 2 B 



V 2 B = /i 



Vx (VxM) + — (VxP) 
at 



(3b) 



Using the Lagrangian (I}, one can show that due to photon-photon scattering induced by the exchange of virtual 
electron-positron pairs, the vacuum will effectively be polarised and magnetised, and the polarisation and magneti- 
sation will be nonlinear to cubic order in the electromagnetic field (see, e.g., Ref. @). These equations can be used 
to analyse the properties of coherent radiation, and it is straightforward to show that the non-linearities vanish in 
the case of co-propagating plane waves. On the other hand, introducing dispersion into the propagation, by the use 
of wave guides, will give non-trivial polarisations and magnetisation, and this has been used as a suggested means of 
detection of photon-photon scattering Q| , and to form light bullets in vacuum |Toj . 

On the other hand, as shown in Ref. one can derive the dispersion relation for a single photon moving on 
a given background of electromagnetic fields, using first principles. The result is (see also Ref. |12| and references 
therein) 



w(Jb,£?,fl)=c|fc|(l-|A|Q| 



(4) 



where 



\Q\ 



so 



c 2 B 2 - (fc • E) 2 - c 2 (fc • B) 2 - 2ck ■ (E x B) 



(5) 



and A = A±, where A+ = 14k and A_ = 8k for the two different polarisation states of the photon. Furthermore, 
fc = fc/fc. The approximation A|Q| 2 <C 1 has been used. The background electric and magnetic fields are denoted by 
E and B, respectively. 

From expressions (@J) and J5J we can derive essentially different dispersion relations, depending on whether or not 
the background is coherent or incoherent. 

To start with, suppose a photon with wave vector fc is moving on a background of photons close to thermal 
equilibrium. The background electromagnetic fields then satisfy (E) = (B) = 0, (EiEj) 
(B 2 )Sij/3, where the angular brackets denote ensemble average. From (JSJ we then obtain 



{E 2 )5 l3 /Z, (BiBj) = 



4 a 
3%' 



where S = Eq((E 2 ) + c 2 (B 2 ))/2 is the energy density of the radiation gas. 

On the other hand, if the photon propagates on a plane wave background, such that E = Ee, and B 
where e is the unit electric vector and k p is the wave vector of the background field, we find that 



(6) 



Ek p x e/c, 



\Q\ 



2 - 2(fc • k p ) - (fc • e) 2 - (fc • (kp x e)f 



(7) 



If fc = fc p , the interaction term \Q\ 2 vanishes identically, consistent with the previous statement that a plane wave 
pulse can not self- interact. On the other hand, if the single photon belongs to a random ensemble of photons, i.e., a 
photon gas close to thermal equilibrium, we have (fc) = and (kikj)Sij /3. The expression (JJJ then reduces to 



\Qp\ 2 — 



(8) 



where & v — e^E^ is the energy density of the background electromagnetic pulse. 
Using the dispersion relation above, different scenarios can be analysed. 
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A. An ensemble of pulses 

For a plane wave pulse moving on a background of photons, randomly distributed in momentum space, the relation 
© holds. It is then straightforward to derive, by using the eikonal representation and the WKBJ approximation [l5| . 
the Schrodinger equation 0] 



c 



V 2 - (fc -Vf E + ^\ck g g E = 0, (9) 



where fco is the linear vacuum wave number of the pulse. Suppose now that we have an ensemble of pulses. We 
decompose the electric field according to E = (E) + E, where (E) = 0. Taking the average of Eq. @, we obtain 



% 



(fco) d , d \ 1/ lt t . \ d 2 (E) 2 ,„..,..„. / fc dE BE 1 1 t d 2 E V n . 



where we have employed Einstein's summation convention. In general, we can subtract Eq. (|10fl from in order to 
obtain a similar equation for E. This would then constitute a general coupled system between the fluctuations and 
the average wave packet. We will not pursue this issue further, however. 

We now write fco = (fco) + fco with (fco) = 0, which implies (fcoifcoj) = (koi)(koj) + (koikoj). If fco is approximately 

isotropically distributed, we then have (koikoj) = |(|fco| 2 )<5y. Furthermore, since |fco| 2 = 1, it follows that < 

(|fc | 2 ) = l-|(fco)| 2 <l. Thus, 

(k 0i k oj ) = (ki)(k 0j ) + | (l - l(feo)| 2 ) Sij. (11) 

In the case of small random fluctuations in the field amplitude, i.e., \E\ <C \(E)\, we neglect the right hand side 
of Eq. (Itjfr. Suppose furthermore that the averaged wave vector, its direction, and its amplitude all change slowly 
over time and space compared to the amplitude (E). Moreover, for a self- interacting pulse ensemble, we see that 
£ — £q(|.E| 2 ) w e |(-B)| 2 . Using the relation iJTTJ, we thus obtain the dynamical equation 



dt (fc > " ' ) x ' 2(fc ) 



2 ^ I , ,>\ 2 



^l(feo)| 2 )vi + |(l-|(fc )| 2 ) Vjj 



{E) + 2XceM | (jE;) |2 (jB) = 0) 



3 (fco) 

(12) 

where = V 2 — V 2 and V 2 | = ((koi)(koj)/\(ko)\ 2 )didj. Equation i|12|) therefore describes the dynamics of the 
non-fluctuating part of an ensemble of pulses, when the nonlinear self-interaction due to photon-photon scattering is 
taken into account. 

B. The kinetic theory of non-linear photons 

In general, from our dispersion relation one can formulate the single-particle dynamics in terms of the Hamiltonian 
ray equations. For a dispersion relation lo = cfc[l — (A/2)|Q| 2 ], where \Q\ 2 is assumed to be independent of fc, we 
have the Hamiltonian ray equations 

(Lc f i - h\Q\A fc, (13a) 



dk V 2 

* = -§H«^' <>») 

where r is the group velocity of the photon, fc is the force on a photon, and the dot denotes time derivative. 
The equation for the collective interaction of photons can then be formulated as 0, 0] 

- + c 1 - -A|Q| 2 fc - - + -Xck^L . _ = o. (14) 



where the distribution function N — N(k,r,t) has been normalised such that the number density is n(r,t) — 
J N(k,r,t) dk. 
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From Eq. (|14f> a hierarchy of fluid equations can be built 01 , of which the two first are the energy conservation 
equation 



Ml 
at 

and the momentum conservation equation 

an 

dt 

respectively. Here 



V • {S g u + q) 



K,9\Q\ 2 



dt 



V • 



u ® n + p 



= -A<fV|Q| 2 , 



hujNdk 



(15a) 



(15b) 



(16) 



is the energy density, q(r,i) — J huwN dk the energy flux, and we have made the split r = u + w, where (w) = 0. 
Furthermore, H = J UkN dk is the momentum density, and P = J w <g> (hk)N dk is the pressure tensor. 



C. The interaction between a radiation gas and an electromagnetic pulse 

From Eqs. © and (JZJ) it is clear that even though the interaction of a pulse with itself will vanish, this is not the 
case if the pulse moves on a background of incoherent photons, i.e., a radiation gas. The problem of an interacting 
pulse and an incoherent radiation gas via photon-photon scattering was investigated in Ref. |l6j |. where a coupled 
system of equations was derived |l5| . 



D. The self-interacting radiation gas 

From JBJ it follows that a kinetic photon gas may interact with itself via a ponderomotive force oc V^ s . Using Eqs. 
© and 114(1 . the kinetic equation for this type of gas is 

dN ( 2 , \ i dN 2 , , d§ dN 

iH +c { l --, xs ) k -^ + r hX ^-i)k- - (17) 

Assuming N(k, r, t) = N (k) + Ni(k) exp[i(K • r - fit)], where N t < N Q and S\r, t)=S Q + S x exp[i(K • r - Qt)}, 
where &\ <C <Sb, we linearise Eqs. I|17fl and (|16fl to obtain the dispersion relation 



i = 2 -xm( k2K ^ dN ^dk. 

3 J n-cK-k 



(18) 



If the background distribution is isotropic, i.e., No(k) — No(k), and introducing spherical coordinates in fc-space, 
we rewrite Eq. p8|l as 

47r, t f 1 a , f°° l4 dN 



l = —\hc da k^—^-dk (19) 

3 7-i P - a J dk 

where (3 = fl/cK and a — k ■ K / K . Performing the angular integration in a, for p < 1, we obtain 

1 = JLxfic [-1 + arctanh {(3) - nip] / fc 4 — - dk. (20) 
3 Jn dk 



We now assume that we have a Gaussian background distribution 



2a 2 



(21) 



where a is the width of the distribution. The dispersion relation (|2U|I is then 



1 = ^A^ [-1 + P arctanh (/3) - iriP) . (22) 
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III. CONCLUSION 



In the present paper we have considered the interaction, due to the quantum electrodynamic photon-photon scat- 
tering, of electromagnetic waves and incoherent photons. For this purpose we have applied the dispersion relation 
derived in Ref. The dynamical equations for an ensemble of pulses has been derived, and a set of equations 

describing a radiation fluid has been presented. An expression showing Landau damping within the self-interacting 
radiation gas has been derived. 
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